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PREFACE 
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research  was  done  under  Project  No.  V37A-32A  in  support  of  the  High 
Reynolds  Number  Tunnel  (HIRT)  Project.  The  author  of  this  report 
was  James  P.  Billingsley,  ARC,  Inc.  The  manuscript  (ARC  control 
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1.0  INTRODUCTION 


It  is  highly  desirable  to  more  accurately  measure  "second-order" 
aerodynamic  effects  resulting,  for  example,  from  Reynolds  number  in¬ 
fluence.  Consequently,  all  wind  tunnel  data  systems  are  required  to 
provide  a  high  degree  of  accuracy.  One  area  that  requires  thorough 
analysis  is  the  magnitude  of  sting  bending  and  oscillations  that  will  exist 
under  high  free-stream  dynamic  pressure  and  impulse  starting  loads 
characteristic  of  some  wind  tunnels  (Ref.  1). 

The  majority  of  aerodynamic  data  available  from  wind  tunnels  is 
obtained  from  captive  (restrained)  model  testing.  In  this  situation,  the 
model  is  supported  at  the  free  end  of  a  beam  (the  sting)  which  is  canti¬ 
levered  forward  from  a  support  sector  into  the  oncoming  airflow  in  the 
test  section.  Aerodynamic  forces  and  moments  are  typically  measured 
with  a  six- component  strain-gage  balance  located  internal  to  the  model. 
The  model  is  attached  to  the  balance,  and  the  balance  is  attached  to  the 
end  of  the  supporting  sting  as  illustrated  in  Fig.  1.  Both  the  sting  and 
the  balance  will  deflect  under  aerodynamic  and  inertial  loads. 


Figure  1.  Simple  sting  deflection  geometry. 
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Many  models  and  stings  will  be  subjected  to  a  rapid  application  of 
load  during  the  tunnel  starting  process  and  will  suffer  large  static  de¬ 
flections  under  certain  important  test  conditions  due  to  the  magnitude 
of  the  load  applied.  Figure  1  illustrates  typical  sting  deflection  geome¬ 
try.  For  certain  test  conditions,  the  sting  tip  angular  and  translational 
deflection  could  be  several  degrees  and  inches,  respectively.  These  de¬ 
flections  are  generally  in  excess  of  sting  deflections  normally  encountered 
in  lower  dynamic  pressure  facilities.  However,  it  is  not'  the  static  de¬ 
flection  alone  which  causes  concern  during  the  portion  of  the  run  when 
data  are  obtained.  The  oscillatory  deflection  of  the  sting,  superimposed 
on  the  static  deflection  which  is  also  present,  can  have  a  major  influence 
on  the  accuracy  of  steady-state  data  (Ref.  2). 

This  investigation  was  undertaken  primarily  to  evaluate  the  struc¬ 
tural  dynamics  of  sting-model  systems  subjected  to  transient  load  con¬ 
ditions.  The  results  are  applicable  to  existing  test  conditions  such  as: 

1.  Blowdown  or  Short -Duration  Facilities: 

The  highly  impulsive  flow  conditions 
create  sting  oscillations  which  can 
influence  data  acquisition  and  accuracy 
and  cause  structural  integrity  loss. 

2.  Facilities  Utilizing  Model  In'jection  Systems: 

The  sting-model  injection  into  the  tunnel 
test  section  under  steady  flow  conditions 
can  create  rather  large  sting  deflections. 

Here,  the  primary  concern  is  for  the 
sting-balance  structural  integrity  rather 
than  data  accuracy. 

3.  Facilities  without  Model  Injection  Systems: 

Tunnel  starting  and  stopping  loads  on  the 
model  and  sting  are  particularly  severe  for 
transonic  and  supersonic  facilities  because 
^of  the  normal  shock  which  traverses  the 
test  section. 

4.  Testing  with  Very  High  Angles  of  Attack: 

High  angle -of -attack  stall  and  buffeting 

of  airplane  models  creates  sting  transient 
deflection  conditions  and  the  associated 
previously  mentioned  problems. 
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5.  Facilities  Utilized  for  Dynamic  Stability 
Testing:  During  dynamic  stability  testing 
for  aerodynamic  damping  derivatives, 
the  basic  sting  oscillatory  motion  can  be 
superimposed  on  the  dynamic  balance 
oscillatory  motion.  This  creates  data 
acquisition  and  accuracy  problems  (Refs. 

3  and  4). 

This  report  discusses  various  sting  materials  and  describes  the 
analytical  techniques  utilized  to  compute  sting  deflections,  spring  con¬ 
stants,  natural  frequencies,  and  transient  motion.  Practical  utilization 
of  these  computational  schemes  is  demonstrated  by  comparing  experi¬ 
mental  model  sting  oscillatory  data  from  a  hypersonic  impulse  tunnel 
with  computed  results. 


2.0  ANALYSIS 


2.1  STING  MATERIALS 

Slender  beam  deflection  is  primarily  caused  by  bending  moments, 
and  the  deflection  is  inversely  proportional  to  El,  where  E  is  Young's 
Modulus  and  I  is  the  cross-sectional  area  moment  of  inertia.  Shear 
deformation  is  usually  less  than  1  percent  of  the  total  deflection,  unless 
the  beam  is  very  short  relative  to  its  diameter  (  L  !=»d  ).  The  allowable 
bending  stress  (ffy)  is  also  a  critical  beam  design  item  which,  for  a 
given  moment,  is  directly  proportional  to  d/I  or  1/d^. 

The  wind  tunnel  test  section  size  limits  the  model  size.  The  sting- 
balance  cross-section  dimensions  are  dictated  by  model  size  limitations 
and  aerodynamic  sting-model  interference  considerations.  Since  the 
sting  eross -section  diameter  is  thus  essentially  specified,  additional 
stiffness  and  lower  stress  levels  cannot  always  be  gained  by  arbitrarily 
increasing  the  diameter  and  moment  of  inertia  (I).  Consequently,  any 
reduction  in  sting  deflection  must  eventually  come  from  inherent 
material  stiffness  parameters,  which  in  this  case  is  obviously  Young's 
Modulus  of  Elasticity  (E). 

Traditionally,  high  strength  alloy  steels  have  been  the  primary 
structural  ingredient  for  wind  tunnel  stings  and  balances.  However,  in 
order  to  combat  excessive  deflection,  other  materials,  which  are  stiff er 
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than  steel,  warrant  serious  consideration.  Table  1  (Refs.  5,  6,  and  7) 
lists  pertinent  characteristics  of  alloy  steels  and  additional  materials 
which  are  potential  candidates  for  sting  construction.  It  is  acknowledged 
that  the  alloy  steels  are  cheaper  and  easier  to  fabricate,  than  any  of  the 
stiffen  materials. 

Table  1  indicates  that  the  sintered  tungsten-titanium  carbide 
alloys  are  about  three  times  stiffen  than  the  steels.  However,  they 
are  somewhat  marginal  so  far  as  the  yield  stress  (o-y)  is  concerned. 

They  also  tend  to  fail  in  a  brittle  manner  when  the  yield  stress  (essen¬ 
tially  the  elastic  limit  stress)  is  exceeded.  Even  so,  they  merit  serious 
attention  as  a  primary  sting  material  because  of  the  threefold  stiffness 
increase  for  cases  in  which  sting  deflections  pose  problems. 

Table  1.  Characteristics  of  Sting  Construction  Materials 


Material 

E  X  10“^ 
Ib/in.^ 

w, 

lb/in.3 

Oy  X  10 
lb/in.2 

Aluminum  Alloys 

10.4 

0.101 

73 

Titanium  Alloys 

16.0 

0.170 

170 

Nickle-Base  Superalloys 

29.0 

0.296 

175 

Steel  Alloys 

27.0  to  30.0 

0.281 

250 

Cobalt-Base  Superalloys 

33.6 

0.320 

116 

Molybdenum  Alloys 

46.0 

0.370 

125 

Tungsten 

59.0 

0.700 

220 

Sintered  Tungsten-Titanium 
Carbide  Alloys 

** 

oo  ‘o*** 
82.2 

0.510 

124**** 

♦ 

Ref.  5 

Ref.  Nominal  Properties 

*S(C5)< 

Ref.  6,  Conservative  Estimate 

^ 

Estimated  from  Weibull's  formula  (Ref.  7) 
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2.2  STING  DEFLECTION,  STRESS,  STRAIN,  AND 
NATURAL  FREQUENCY 


Because  of  model  space  limitations  and  aerodynamic  sting-model 
interference  effects,  sting  configurations  normally  consist  of  one  or 
more  linearly  tapered  sections  combined  with  one  or  more  untapered 
sections.  The  balance  strain-gage  lead  wires  are  generally  brought 
out  through  a  hole  down  the  center  of  the  sting.  Thus,  cross-section 
dimensions  will  usually  not  be  uniform  throughout  the  length  of  the  sting. 
These  numerous  cross-section  discontinuities  make  closed  form  inte¬ 
gration  of  the  basic  slender  beam  elastic  curve  equation  (Eq.  (1))  im¬ 
practical  for  realistic  stings.  The  basic  slender  beam  elastic  curve 
equation  is 

d^z  _  M  • 

-  El  (1) 

where  M  is  the  total  bending  moment  at  x,  and  E  and  I  are  usually  func¬ 
tions  of  X  also.  Consequently,  numerical  integration  of  Eq.  (1)  is  re¬ 
quired  for  practical  sting  configurations.  Because  stings  are  normally 
cantilevered  from  some  rigid  structure,  they  are  ideally  suited  for 
analysis  by  the  first  and  second  area  moment  propositions  (Ref.  8).  The 
area  moment  analysis  was  incorporated  in  a  computer  program  which 
was  written  for  the  cantilever  beam  schematically  illustrated  in  Fig.  2. 
This  beam  has  four  different  cross-section  discontinuities.  Each  section 
can  either  be  a  constant  diameter  cylinder  or  a  linearly  tapered  cone 


El  E2  E3  E4  E5 


Notes;  R(N1)— R(NEND)  =  Radius  at  Section  Discontinuities 

N1---NEND  =  Station  Numbers  at  Section  Discontinuities 
Wl— W5  =  Weight/Volume  of  Material  in  Each  Section 
El— E5  =  Young's  Modulus  of  Material  in  Each  Section 
RHl— RH5  =  Radius  of  Hole  through  Section 

Figure  2.  Beam  geometry  and  basic  input  for  STING-1  program. 
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frustum.  Each  section  may  have  a  hole  of  arbitrary  constant  diameter 
along  the  centerline.  The  material  of  each  section  is  also  arbitrary 
(i.  e.  steel,  aluminum  alloys,  carbide  alloys,  etc.).  The  cross- 
section  shape  must  be  a  circle  or  concentric  circles,  although  this  is 
not  a  requirement  of  the  basic  beam  analysis. 


This  computer  program,  STING-1,  can  perform  the  following 
calculations  after  the  basic  geometric  and  material  properties  are 
provided. 


1.  Compute  deflections  (angiTar  and  translational)  along  the  beam 
length  caused  by  the  weight  of  the  beam.  This  information  is  the  static 
deflection  curve  of  the  beam  alone  and  is  necessary  for  the  computation 
of  the  first  natural  bending  frequency  by  Rayleigh's  Method  (Eq.  (4)), 
Maximum  bending  stress  and  strain  (along  the  beam  length)  caused  by 
beam  weight  are  also  computed  by  the  following  standard  relations, 
respectively: 


Me 

I 


(2) 


(3) 


2.  Deflection,  stresses,  and  strain  caused  by  concentrated  loads 
or  moments  at  any  station  along  the  beam  may  be  computed.  In  addition 
to  obvious  design  applications,  the  deflections  caused  by  concentrated 
loads  are  necessary  to  compute  the  beam  spring  constants  or  influence 
coefficients  (see  Section  2.  3).  Figure  3  shows  bending  stresses  for 
a  sting  with  a  concentrated  load  of  24,  000  lb  at  one  end  computed  by  Code 
STING- 1.  Table  2  lists  deflections  for  steel,  carbide,  and  a  composite 
carbide-steel  (C-S)  sting  computed  by  STING- 1  for  these  same 
conditions. 


3.  The  first  natural  frequency  may  be  computed  by  Rayleigh's  Meth¬ 
od  (Ref.  9).  The  effect  of  concentrated  masses  (the  model  for  instance) 
can  also  be  included.  The  expression  for  the  sting  natural  frequency  is 


2 


tui 


f  dx 

o  El 


L  2 

J  z'^dm  .+ 
o 


n 

s 


i=l 


g 


(4) 


where  M  is  the  total  bending  moment  (including  beam  weight  and  con¬ 
centrated  weights)  at  station  x,  z  is  the  static  deflection  at  station  x,  z^ 
is  the  static  deflection  at  the  concentrated  weight  locations  (xj^),  and  dm 
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is  the  differential  mass  of  the  beam  at  station  x.  When  the  static  de¬ 
flection  curve  is  used  in  Eq.  (4),  the  computed  frequency  will  be  slightly 
higher  than  if  the  true  dynamic  deflection  curve  is  used.  This  occurs 
because  the  small  deviation  of  the  static  deflection  from  the  first  mode 
dynamic  deflection  imposes  some  constraint  which  makes  the  beam  more 
rigid  and  its  frequency  higher.  This  is  illustrated  in  Fig,  4  where 
frequencies  for  the  carbide-steel  (C-S)  sting  are  given;  both  for  the 
Rayleigh  method  with  static  deflection  and  the  analysis  of  Section  2. 4 
where  the  true  dynamic  mode  deflection  shape  is  computed,  A  compari¬ 
son  shows  that  the  Rayleigh  method  frequency  is  approximately  6  per¬ 
cent  greater  than  the  true  djmamic  mode  frequency.  Figure  4  also  illus¬ 
trates  the  influence  of  material  and  model  weight  on  sting  natural 
frequency.  The  natural  frequency  of  mono-material  beams  (without 
concentrated  weights)  scales  according  to: 


Using  material  properties  from  Table  1, 

^  C  arbi  de 
‘^Steel 


All  Dimensions  Are  in  I  nches 


Station,  x,  in. 


Figure  3.  Maximum  static  bending  stress  in  the  sting. 
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Table  2.  HIRT  Sting  Deflections  and  Natural  Frequencies 


Material 

End  Rotation, 
deg 

End  Deflection, 
in . 

Frequency, 
(Fig.  4,  Wm 

cps 
=  0) 

Steel 

7.37 

8.30 

17.0 

Carbide 

2.47 

2.74 

21.5 

Composite 
(Steel  and 
Carbide) 

2.88 

3.81 

16.0 

To  compute  the  above  information,  the  following  material  proper¬ 
ties  were  used: 

Material  Ib/in .  Ib/in . 

Steel  30(10®)  0.281 

Carbide  90(10®)  0.525 


Model  Weight,  W|vi,  lb 

Figure  4.  Lowest  natural  frequency  of  the  example  sting  for 
various  model  weights  and  sting  materials. 
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The  ratio  from  the  more  exact  Eq.  (4)  (Wjyj  =  0)  is: 

=  1.26 

‘^Steel  ^'^•0 

4.  Certain  structural  dynamic  analyses  require  the  beam  to  be 
considered  as  a  series  of  lumped  masses  connected  by  springs.  The 
selection  of  a  particular  beam  portion  or  section  to  form  a  concentrated 
mass  is  a  decision  requiring  some  discretion  for  best  results.  Once 
the  partitioning  selection  is  made,  the  STING- 1  program  will  compute 
the  value  of  the  lumped  mass  (mj^),  its  mass  moment  of  inertia  (lyj^)j 
and  its  centroidal  position  (x^).  This  centroidal  position  is  the  point 
at  which  concentrated  unit  loads  and  moments  will  normally  be  applied 
to  compute  the  beam  spring  constants  or  influence  coefficients. 

2.3  INFLUENCE  COEFFICIENTS  (SPRING  CONSTANTS) 

The  structural  dynamic  analysis  described  in  the  next  two  sections 
represents  the  sting-balance -model  system  as  a  series  of  "lumped 
masses"  or  concentrated  masses.  These  masses  are  considered  to  be 
connected  together  by  flat  weightless  springs.  A  load  applied  to  one 
mass  will  influence  the  deflection  of  the  other  masses  because  of  the 
spring  connections.  Consequently,  the  associated  spring  constants 
have  become  known  as  structural  influence  coefficients. 

References  10  (pp.  17-22),  n  (pp,  19-27),  and  12  contain  excellent 
discussions  on  structural  influence  coefficients  which  basically  can  be 
classified  in  one  of  two  categories  as: 

1.  Flexibility  influence  coefficients,  defined  by: 

N 

q;  =  2  c..  Q.  (6a) 

*  j=i  J 

where 

N  =  number  of  generalized  forces  applied 

qj  =  deflection  (translational  or  rotational) 
at  point  i 

Qj  =  Forces  (or  moments)  applied  at  point  j 
Cjj  =  flexibility  influence  coefficients 
=  deflection/unit  load 
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Equation  (6a)  can  be  written  in  matrix  form  as  follows; 


^1 

^2 


^12  ■  ■  ■  ^IN 

f-H 

O' 

Q2 

C]vj2 

Qn 

(6b) 


Basically,  Cj_j,  corresponding  to  the  location  of  the  lumped  masses 
of  the  beam,  can  be  computed  from  the  deflection  (translational  and 
rotational)  caused  by  unit  concentrated  loads  or  moments  applied  at 
the  various  lumped  mass  eg  locations  {x]) . 


For  computational  purposes,  it  is  important  to  note  the  four  differ¬ 
ent  types  of  elements  which  make  up  the  matrix  (see  Ref.  10,  p.  22). 
The  Cy  matrix  can  be  partitioned  as  follows: 


CTF.. 

fcTM.. 

L 

n 

L 

CAF.. 

CAM.. 

n 

i]_ 

N  =  2n 


(7) 


where 


CTFjj  =  translational  deflection  at  i  caused  by  unit 
force  at  j,  ft/lb 

CTM-  =  translational  deflection  at  i  caused  by  unit 
moment  at  j,  ft/(ft-lb)  =  1/lb 

CAFij  =  angular  deflection  at  i  caused  by  unit  force 
at  j,  rad/lb  or  1/lb 

CAM'-  =  angular  deflection  at  i  caused  by  unit  moment 
at  j,  rad/(ft-lb)  or  l/(ft-lb) 

This  allows  C^j  to  be  computed  by  systematic  application  of  the 
STING- 1  program  to  each  of  the  concentrated  mass  positions. 


2.  Stiffness  Influence  Coefficients,  given  by 


(8a) 
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where 


qj  =  deflection  (translational  or  rotational)  at 
point  j 

Qj  =  force  (or  moment)  applied  at  point  i 
K|j  =  stiffness  influence  coefficients 
=  load/deflection 


In  matrix  form. 

Eq. 

(8a)  is 

written 

as 

Qi 

’Kii 

K12  . 

1 

1— 1 

Q2 

= 

K21 

K22 

K2N 

Qn 

Kn2 

^NN. 

'll 

^12 

% 


(8b) 


Note  that  is  completely  analogous  to  the  spring  constant,  K  =  lb/ in. 
employed  in  one-dimensional  spring  analysis. 


The  Ky  coefficients  are  difficult  to  compute  directly  from  beam 
theory.  Consequently,  they  are  generally  computed  by  inverting  the 

This  is  the  manner  K;[j  was  obtained  in  the  present  analysis 


Cj^j  matrix. 


(9) 


Another  important  property  of  the  Cy  and  matrices  is  that  both  are 
symmetric.  The  Clj  matrix  as  computed  by  STING- 1  was  not  symmet¬ 
rical  after  the  third  or  fourth  decimal  places  (round  off  error,  approxi¬ 
mation,  etc. ).  Thus  Cj^j  had  to  be  made  symmetric  so  that  it  would 
properly  invert  and  yield  the  following  multiplication  check: 


fii][s]=[>]  (10) 

where  [I]  is  the  identity  matrix 

It  is  important  to  note  that  the  coefficients  can  be  determined 
experimentally  (Ref.  12).  In  fact,  this  is  preferable  to  computing  them 
analytically  and  may  be  absolutely  necessary  for  complex  structures. 
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In  any  case,  some  experimental  Cij  data  are  valuable  for  verifying  and 
complementing  theoretical  results.  The  VKF  balance  calibration  labo¬ 
ratories  currently  measure  deflection  under  loads  applied  at  the  model 
attachment  point  so  that  static  sting  rotation  from  model  airloads  may 
be  accounted  for  in  data  reduction.  This  type  of  experimental  informa¬ 
tion  was  essential  to  the  application  which  will  be  discussed  in  Section  3.2. 


2.4  STING  NATURAL  FREQUENCIES  AND  MODE  SHAPES 


The  first  natural  bending  frequency  may  generally  be  computed  with 
sufficient  engineering  accuracy  by  the  Rayleigh  method  as  indicated  in 
Section  2.2.  For  design  purposes,  primarily  to  avoid  resonance  condi¬ 
tions,  the  higher  natural  frequencies  are  occasionally  required.  Refer¬ 
ence  9  (pp.  169-170)  outlines  a  procedure  by  which  the  influence  coeffi¬ 
cient  (CTFj[j)  in  conjunction  with  the  lumped  mass  concept  can  be  uti¬ 
lized  to  compute  the  first  and  higher  natural  steady-state  bending  fre¬ 
quencies  and  mode  shapes.  From  Section  2.3,  the  CTFj^j  coefficients 
are  defined  as  translational  normal  deflections  at  x^  caused  by  a  unit 
force  applied  at  xj.  The  final  matrix  formulation  is: 


CTF^2™2 

•  •  CTF. 

In  n 

CTF2imi 

CTF22ni2 

••  CTF2,m^ 

CTF^jm^ 

CTFn2ni2 

•  •  CTF^^m^ 

(11) 


n 

where  A=  1/w  .  This  is  the  standard  matrix  eigenvalue  problem  where 
A  is  the  eigenvalue  and  the  associated  eigenvectors  (zj^)  define  that  par¬ 
ticular  normalized  mode  shape.  Once  the  [CTFymj^]  matrix  is  formu¬ 
lated,  the  eigenvalues  (frequencies)  and  associated  eigenvectors  (nor¬ 
malized  mode  shapes)  may  be  computed  by  the  MATVEC  program  (Ref. 
13).  Concentrated  weights  such  as  the  model  may  also  be  taken  into 
account  by  proper  modification  of  the  appropriate  m^  (see  Figs.  4  and  5 
for  MATVEC  results  for  the  C-S  sting).  The  mode  shapes  illus¬ 
trated  in  Fig.  5  are  typical  for  cantilever  beams. 


2.5  LAGRANGIAN  FORMULATION  OF  STING-MODEL 
MOTION  EQUATIONS 

Sting-balance -model  systems  tested  in  blowdown  or  short-duration 
wind  tunnels  are  exposed  to  transient  aerodynamic  loads  which  can  act 
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on  both  model  and  sting.  Since  useful  data  taking  time  is  limited, 
modern  blowdown  facilities  are  designed  to  vary  the  sting  pitch  angle 
(np)  during  the  run  so  that  the  quantity  of  data  obtained  is  increased 
Continuous -flow  tunnels  may  also  have  this  capability  for  similar  reasons. 
In  addition  to  sector  transient  rotation  capabilities,  some  modern  faci¬ 
lities  are  capable  of  translating  the  sting-model  from  a  location  outside 
the  test  section  into  the  flow  field  of  the  test  section  via  a  so-called 
"injection"  system. 


Figure  5.  Basic  bending  rnode  shapes  for  the  example 
carbide-steel  sting. 

Consequently,  any  realistic  analysis  of  sting  structural  dynamics 
must  account  for  elastic  deflection  of  the  system  caused  by: 

1.  Transient  air  loads  on  both  model  and  sting. 

2.  Transient  loads  associated  with  rigid  body  motion 
of  the  system. 

To  investigate  the  general  transient  planar  motion  of  a  sting-model 
system,  it  has  been  represented  by  a  series  of  concentrated  masses 
as  illustrated  in  Fig,  6a.  These  concentrated  masses  are  connected 
by  springs  so  that  the  motion  of  one  mass  influences  the  motion  of  the 
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other  masses.  Because  of  the  rigid  body  motion  of  the  sting,  it  is 
necessary  to  employ  three  basic  axis  systems  which  are  illustrated  in 
Fig.  6: 

Xj,  Zj  Inertial  axis,  Xj  is  parallel  to  V^. 

y  y 

^IP  Th  is  axis  is  fixed  to  the  nndeflected  sting 

and  moves  with  it  as  shown. 

Xju  ,  Zjjj  This  axis  is  fixed  to  the  ith  concentrated 
^  *  mass  and  moves  with  it  as  shown. 

Lagrange's  formulation  of  dynamics  (Refs.  14  and  15)  is  employed 
to  derive  the  equations  of  motion  for  the  mass.  Lagrange's  method 
provides  a  methodical  and  systematic  approach  which  is  especially 
valuable  in  the  present  application  involving  two  moving  non-inertial 
coordinate  systems. 

The  present  formulation  is  an  extension  to  the  example  of  the 
cantilever  beam  and  single  concentrated  mass  given  in  Ref.  11  (pp.  40- 
41).  The  present  analysis  is  also  similar,  in  some  respects,  to  the 
classical  example  of  airfoil  flutter  given  in  Ref.  11  (pp.  210-212)  and 
Ref.  10  (pp.  532-536)  where  Lagrange's  formiilation  is  employed. 


a.  Planar  deflection  geometry  of  the  sting-model  concentrated  mass  representation 

Figure  6.  Sting  coordinate  systems. 
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2.5.1  General  Lagrange  Equation 


The  Lagrangian  formulation  of  dynamics  provides  a  systematic 
approach  that  gives  directly  the  equations  of  motion  in  whatever  coor¬ 
dinate  system  may  be  convenient  (Ref.  14).  Application  of  Lagrange's 
procedure  is  remarkably  straightforward  even  for  relatively  complex 
systems  involving  moving  non-inertial  coordinate  systems.  Lagrange's 
basic  equation  is 


d  /dT\  _ 
dt\dqj  dq^ 


au 

aqk 


Qk 


Inertial  forces 
(or  moments) 


Conserva¬ 
tive  forces 
(or  moments) 


N  on- conservative 
forces  (or  moments) 


(12) 


where; 

k  =  1— N 
N  =  2n 

N  =  number  of  degrees  of  freedom  or  independent 
coordinates 

n  =  number  of  lumped  masses  in  the  system 
representation 

=  generalized  independent  coordinate.  In  the  present 
analysis,  q^.  will  either  be  angular  rotation  or  a 

linear  translation  distance,  both  of  which  are  easily 
visualized  physical  quantities.  The  motion  direc¬ 
tion  specified  by  Eq.  (12)  is  in  the  exact  direction 
associated  with  q^  (Ref.  15,  pp.  157-158) 

T  =  total  kinetic  energy  of  system  and  must  be  mea¬ 
sured  relative  to  an  inertial  axis  system 
(Ref.  14,  p.33) 

U  =  potential  energy  of  the  system.  In  the  present 
application,  U  is  the  work  done  by  the  sting 
acting  as  an  elastic  spring. 

Qk  =  non-conservative  generalized,  applied  force 

(or  moment)  acting. on  the  ifLmass,  (m^)  corre¬ 
sponding  to  the  ktb  generalized  coordinate.  In 
this  application,  Qk  is  primarily  of  aerodynamic 
origin. 

t  =  time,  sec 


The  appropriate  formulations  of  T,  U,  and  are  given  in  the 
following  sections. 
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2.5.2  The  Kinetic  Energy  and  Its  Derivatives 


In  applying  Eq.  (12)  to  the  sting-model  system  of  concentrated 
masses  as  represented  in  Fig.  6,  a  very  crucial  item  is  the  evaluation 
of  the  total  kinetic  energy  and  its  associated  derivatives.  These  yield 
the  inertial  forces  (or  moments)  of  the  system.  The  total  kinetic 
energy  for  this  system  is 

n 

(n  =  number  of  masses)  (13  a) 


T  =  X  T; 
i=l  ^ 


Ti  =  2/  V 
m  • 


abs 


(13b) 


Tj^  is  the  kinetic  energy  of  the  mass  (mj)  and  is  the  absolute 

velocity  of  an  element  of  mass  (dm),  which  is  located  as  shown  in  Figs. 
6b  and  c  with  respect  to  the  inertial  axis  (Xi,  Zj),  The  coordinates 
of  the  particle  (dm)  with  respect  to  the  Xj,  Zj  origin  are 


X]  =  +  Xjj  cos  ctp  +  Zjj^ 


+  Xjjj  cos  (a„  +  0;)  +  Zjjy  sin  (a„  +  d-) 


III. 


\  =  “p  +  “p 


—  Xjjj  sm 


in  (a„  +  0;)  +  Zjjj  cos  (a„  +  $■) 


(14  a) 
(14b) 


Local 
Deflected 
Sting  Axis 


b.  Coordinates  of  mass  element  with  respect  to  moving  axis  systems 

Figure  6.  Continued. 
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c.  Coordinates  of  mass  element  with  respect  to  an  inertial  axis  system 

Figure  6.  Continued. 


Differentiating  Eq.  (14)  with  respect  to  time  yields  the  absolute 
velocity  components  of  the  particle  (dm).  Note  that  and  zjjj.  =  0 

since  the  axis  is  fixed  relative  to  mass  (mp  and  rotates  with 

it  through  the  angle  +  ©i).  The  origin  of  the  Xjjj. ,  Zjjj.  axis  is  not 
necessarily  the  center  of  gravity  of  the  mass  (m^^)  (Fig.  6cL  Note  also 
that  xjj.  is  a  constant  distance  which  denotes  the  location  of  the 
axis  origin  along  the  sting  length.  Strictly  speaking,  xjj^  does  vary  slight¬ 
ly  as  the  sting  deflects  under  load,  but  this  very  small  change  will  be  ne¬ 
glected,  The  final  relation  for  the  total  kinetic  energy  is 


T 


■2 


i=l 


+  2q.(xQsina 


2ko  Xjip  sinap- 


COS 


+  Xj  6 


2 

P 


+  cosap  -  xjip)  +  q?  +  2q.(k^apCosap 
+  q?«p]mi  -  [x„  sin  (ap  +  q^^^)  +  cos  {a^  +  qj^J 


+  (qi  "  xjflp)  cos(q.^J  +  q.  dp  sin 

+  [k„  cos  (ap  +  q^^J  -  z^sin  (ap  +  q^^^)  -  {q.  -  x.  dp)  sin  (q^^^) 

+  qi  «P  cos  (q.^^)](dp  +  qi+n)(zimi)  +  '/2(dp  +  (15) 
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where 

=  ^iij  =  ^ 

^i+n  = 

'‘i  =  ^iij 

r  dm  =  m. 

Jm .  1 

1 

fm.  ’^III.  =  X.  mj 

/ra.  ^Ul.  "’i 

/m.  +  ^ni)  =  ^y; 


(see  Fig.  6d) 

(16a) 

(see  Fig.  6d) 

(16b) 

(see  Fig.  6d) 

(16c) 

(17a) 

(see  Fig.  6c) 

(17b) 

(see  Fig.  6c) 

(17c) 

(17d) 

Note  that  Xq  and  Zq  are  the  translational  velocities  of  the  model  support 
system  moving  as  a  rigid  body.  This  motion  will  occur,  for  instance, 
when  a  model  is  injected  into  a  wind  tunnel  on  a  system  in  which  the 
model,  sting,  and  support  translate  in  unison. 


d.  Final  coordinates  of  deflected  sting  with  respect  to  undeflected  sting 

Figure  6.  Concluded. 


24 


AEDC-TR-76-41 


Straightforward  evaluation  of  the  kinetic  energy  derivatives  yields 
the  following  results: 

Force  Inertial  Terms  (k  =  i) 

d/5T\  <9T  ...  ..  •  2  ••  ^ 

lli~)  ~  ^  ~  +  x^sinap)mj 

\^q./  ^1 

-  [cos  (q.^^)Cip  +  ki+r)  -  (18a) 

-  [sin(qi^J(ap  +  -  cos  (q.^JCSp  +  qi+J mj) 

Moment  Inertial  Terms  (k  =  i  +  n) 


(9t 


=  (qi+n  +  «p)iyi  -  t(qi  + 
+  +  2  qi  ip  +  q;  ip 


«p  “  <ii  (qi+n) 

^  X.  i^)  sin  (q.^^)](V.  m.) 


-  [(qi  +  -  ^i«p  "  qi«p)  (qi+n)  (^o  +24.^ 

+  qi  «p  -  ^i  ^)  (qi+n)](^i  "'i) 


2.5.3  The  Potential  Energy  and  Its  Derivatives 


In  Eq.  (12),  the  partial  derivative  (au/9qj^)  represents  that  part  of 
the  generalized  forces  that  can  be  derived  from  a  potential  function  (U). 
In  the  present  system,  these  are  the  stabilizing  restoring  forces  created 
by  the  elastic  deformation  of  the  sting.  The  sting's  capacity  to  do  work 
is  equal  to  its  stored  potential  strain  energy,  which,  according  to  Ref. 

11  (p.  22),  can  be  written  as  follows  in  terms  of  the  stiffness  influence 
coefficients : 


U 


■  i 

=  ^4[Kjj  qj  +  K21  qg  qi  +  -- -  -  Kfjj  qjj  q^ 

+  K;^2  qi  q2  +  ^22  qi  +  "  ■  ■  ^n2  qN  q2 


(19) 


+ 


+  Kin  qi  qN  ^2N  q2  qN  +  ^nn 
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It  is  shown  in  Ref.  11  that 


-  Kji  +  Kj2  42  +  -  -  -  Kjjyj  4]^ 

=  Qi 


=  Force  (or  moment)  at  point  i  in  the 
direction  of  q.  (see  Eq.  (7a)) 


(20) 


In  the  present  formulation, 


dq- 


2n 


(21a) 


Sting  elastic  restoring  force  at  point  x. 


dU  _  2^1 
j=l 

=  Sting  elastic  restoring  moment  at  point  Xj 


(21b) 


2.5.4  The  Non-Conservative  Generalized  Forces 

Equation  (12)  provides  for  the  arbitary  representation  of  time  varying 
generalized  forces  (Qj^)  that  may  drive  and/or  damp  the  system  motion. 

Qk  is  a  non- conservative  force  or  moment  which  essentially  means  that 
it  is  not  derived  from  a  potential  function  such  as  potential  energy  (Sec¬ 
tion  2.  5.  3).  For  example,  dissipative  effects  such  as  sting-model  aero¬ 
dynamic  damping  and  sting  structural  damping  are  clearly  in  this  cate¬ 
gory. 

2.5.4. 1  Unsteady  Aerodynamics 

In  the  present  application,  the  primary  driving  force  (Qk)  is  creat¬ 
ed  by  unsteady  aerodynamic  phenomena.  The  resulting  aerodynamic 
loads  may  be  obtained  by 

1.  Theoretical, analysis  -  Refs.  10  and  11  discuss  and 
critique  the  available  theoretical  efforts.  Basic 
trends  are  reasonably  well  defined  for  ideal  config¬ 
urations. 
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2.  Experimental  measurements  -  Good  experimental 
data  are  highly  desirable  and  should  be  utilized 
when  available  (see  Section  3.2). 

3.  Approximate  analysis  utilizing  both  theoretical  and 
experimental  information  -  Practical  requirements 
for  unsteady  aerodynamic  simulation  of  realistic 
configurations  may  dictate  such  an  analysis. 

The  following  quasi-steady  approximations  are  used  in  the  present 
formulation.  For  k  =  i,  the  aerodynamic  normal  force  is 


Qi  =  N.p  =  ^  S;  C. 


T. 

1 


(22) 


where  q^  may  vary  with  time  and  C]^,p.  is  composed  of  both  a  "static" 
and  "dynamic"  contribution  as  follows; 


1  ;  11  a 


Static  Force 
Coefficient 


Dynamic  Force 
Coefficient 


(23) 


If  (o'Tj^)  is  linear,  then 

Cn  ^  ~  (24) 

1  i  “i  * 

If  nonlinear,  then  C]^ .  is  curve  fitted  to  a  suitable  function  of  o-rp..  In 
a  similar  manner  for  k  =  i  +  n,  the  aerodynamic  moments  are  ^ 


where 


then. 


Qi+n  =  -  (zj  -  z.)] 


(25) 


M 


T.  ~  4x)  °i 


S.  d,  c. 


aerodynamic  pitching  moment  of  mass 
(mj)  at  point  (xj) 


(26) 


T. 


a-  \  2V  /  m.,  \2V  / 

oo/  Mj\oo/ 


(27) 


Slali  c 
Moment 


Dynamic  Damping  Moment 
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Also,  if  J  is  linear,  then 


~  “T. 

11  aj  1 


(28) 


Otherwise,  it  is  curve  fitted  as  a  function  of  Now  instan¬ 

taneous  angle  of  attack  at  station  xi  along  the  sting,  is: 


Also: 


'T. 

1 


«p  +  q 


+ 


1  P 


“  li 


i-l-n  y  cos  a 


Angular 
Rotation  . 


Plunging  Translation 


(Xi^  -  qj)V^ 


^  i+n  y  CO s  a 


(xia  -  4j)a  sin  a 


V  cos  ct. 


cos 


V  z  -  V 

CO  O  oo  O 


The  rotational  part  of  o-Xj^  is  which  is 


*Ir.  +  ^i-t-n 


=  Gp  +  0i 


AF^  is  the  instantaneous  axial  force  on  the  model  (mj^),  which  is 

=  900^1  Cat 


(29) 


(30) 


(31) 


(32) 


The  axial  force  contribution  is  a  stabilizing  influence  (see  Fig.  7)  with 
respect  to  the  beam  since  its  moment  apparently  acts  to  reduce  deflec¬ 
tion,  at  least  for  the  first  mode  shape.  In  most  cases,  the  axial-f.orce 
moment  is  largest  along  the  aft  and  stiffen  portion  of  the  sting  so  that 
the  stabilizing  influence  is  minimized. 


The  above  formulation  for  (Eq.  (23))  and  (Eq.  (27))  is 

called  "quasi-steady"  because  it  implies  instantaneous  response  to  ar£. 
This  representation  has  been  used  in  air  craft -missile  trajectory 
computer  programs  (Refs.  16  and  17,  for  instance)  with  reasonably 
satisfactory  results.  It  has  also  been  used  in  Ref.  2  for  sting-model 
motion  simulation.  Consequently,  its  use  is  recommended  in  the 
present  application  unless  more  precise  theoretical  and/or  experimen¬ 
tal  information  is  available.  Certain  configurations  and  test  conditions 
may  require  a  more  sophisticated  aerodynamic  force  and  moment 
representation. 
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Inclusion  of  the  angular  rate  (ur^.  and  qp.)  dependent  aerodynamic 
forces  and  moments  has  a  negligible  effect  on  motion  amplitude  and 
decay  for  the  relatively  stiff  stings  commonly  employed  in  standard 
static  stability  and  control  testing.  However,  attention  must  be  focused 
on  their  influence  in  dynamic  stability  testing  when  a  highly  sensitive 
dynamic  balance  is  utilized. 


Sting  at  xi  =  0 

Axial  Force  Moment  =  AFi  [xj0i  -  (zi  -  zj)  ] 


Figure  7.  Axial  force  bending  moment  schematic. 

2.5.4.2  Structural  Damping 

Observation  of  experimental  sting  motion  data  obtained  from  the 
VKF  Tunnel  F  (Section  3. 2)  revealed  that  the  motion  was  highly  damped. 
This  timewise  decrease  in  amplitude  was  much  more  than  could  be 
accounted  for  by  aerodynamic  damping  alone.  An  investigation  revealed 
that  the  decrease  in  amplitude  was  essentially  linear  with  increasing 
time.  This  manner  of  amplitude  decay  is  normally  associated  with 
Coulomb  or  friction  damping  (Ref.  9)  and  the  single  degree  of  freedom 
differential  equation  is  of  the  form: 


mx  + 


(33) 
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F  is  a  constant  friction  force  and  always  acts  to  oppose  the  motion. 
Consequently,  it  was  suspected  that  some  form  of  structural  damping 
capable  of  dissipating  a  large  amount  of  energy  was  prevalent  in  the 
VKF  Tunnel  F  sting. 

A  brief  literature  survey  revealed  that  there  are  basically  two 
types  of  structural  damping: 

1.  Internal  damping  of  a  solid  homogeneous  structural 
member  caused  by  internal  grain  structural,  plastic 
strain,  magnetic-electrical,  or  thermal  effects 
(see  Ref.  18  for  a  comprehensive  survey),  and 

2.  Damping  of  joints,  connections,  or  interfaces  of 
the  members  of  a  total  structure  caused  by  inter - 
facial  Coulomb  friction  and  slip  (see  Refs.  19  and 
20). 


Reference  20  reports  that  under  optimum  circumstances,  the 
energy  dissipated  in  joint  slip  damping  may  be  large,  not  only  compared 
with  internal  damping,  but  also  compared  with  the  total  strain  energy 
of  the  structure.  Since  the  VKF  tunnel-sting  contained  several  promi¬ 
nent  joints  and  its  motion  was  highly  damped  with  essentially  a  linear 
decay,  it  is  concluded  that  damping  is  primarily  caused  by  interfacial 
Coulomb  friction. 


Consequently,  it  is  postulated  that  the  sting  damping  coiild  be 
modeled  in  a  rather  gross  manner  by  the  following  relations  which  are 
somewhat  analogous  to  the  Coulomb  friction  force  in  Eq.  (33).  The 
damping  force  at  xi  is 


QDFj  =  -FCOFj 


2n 

2 

i=i 


(34a) 


The  damping  moment  at  is: 


where 


and 


QDMj  =  -FCOMj 


i+n,]  1] 


2n 

2 

J=i 


(34b) 


2n 


2 

j=i 


‘‘i+n.i  ’j 
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are  the  sting  elastic  restoring  force  and  moment,  respectively,  at  Xj^ 
(Section  2.  5.  3).  FCOFj|^  and  FCOMj^  are  arbitrary  constants  which  are 
utilized  to  adjust  the  fractional  magnitudes  of  the  sting  elastic  restoring 
forces  and  moments,  respectively,  which  are  resisting  the  sting  motion. 
They  are  determined  by  numerical  experimentation  and  comparison 
with  experimental  sting  motion  damping  data  (see  Section  3.  2). 

Formulation  of  a  more  elaborate  mathematical  model  of  sting 
structural  damping  than  Eq.  (34)  will  require  both  analytical  and  experi¬ 
mental  results  for  typical  wind  tunnel  sting  joints. 

2.5.5  Specific  Motion  Equations 

The  specific  motion  equations  of  each  individual  concentrated  mass 
(mj^)  of  the  sting-model  system  for  each  of  its  two  degrees  of  freedom 
corresponding  to  the  general  formulation  (Eq.  (12))  can  now  be  written 
by  combining  appropriate  results  from  Sections  2.5.2,  2.5.3,  and 

2.  5.4. 

For  k  =  i,  Eqs.  (18a),  (21a),  (22),  and  (34a)  are  substituted  in  Eq. 
(12)  to  yield  the  force  relation  (in  the  translational  Zj^  direction)  for  mass 
(mi).  It  is  convenient  for  future  purposes  to  write  this  equation  as 
follows: 


where: 


+  ""12.  =  (RHSF)j 

mil,  =  mi 


(35a) 
(36  a) 


mi2.  =  -  [cos  (qj^^){x.  mj)  +  sin  (qj^^)(z.  nii)]  (36b) 

(RHSF)i  =  (xjCp  +  qja^  -  z^cosap  -  i„sinap)ini 

+  [cos(qj^J(zj  mj)  -  sin(qi^^)(xi  m^Kqj^^  +  ap)^  {36c) 

2n 

-  Kij  qj  +  N-p^  +  QDFj 


For  k  =  i  +  n,  proper  combination  of  Eqs.  (18b),  (21b),  (25)  and 
(34b)  yields  the  moment  equation  (in  the  rotational  direction)  for  mass 
(mi).  This  is: 
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""21.  +  ‘^22/qi+n  +  =  (RHSM)i 


(35b) 


where 


^21.  —  ^12-  (36d) 

m22.  =  ly.  (36e) 


(RHSM).  =  [(z^  -  X.  a'p  -  qi  cos  (qj^^)  +  +  2qj 

+  qi  «p  -  «p)  sin  (qi+n)](xi  mj)  +  [{z^  -  x-  flp 

-  +  2qi  Op  +  qi  Hp  -  X.  a^)  cos  (qi+JKzi  m.) 

-  AFj[x.(qj^J  -  (qj  -  qj)]  +  M^,  +  QDMj 


qi  Op)  sin  (qi+n)] 

2n 

jJl  ^i+n,j 


(36f) 


Equation  (35)  is  comprised  of  the  basic  equations  which  are  both 
statically  and  dynamically  coupled  together  in  the  general  case.  Nrj-^ 
and  Mq-.  are  still  in  arbitrary  form.  If  N-p.  and  Mt.,  as  represented 
byEqs.^(22),  (23),  (26),  (27),  (29),  (30),  and  (31),  ^e  substituted 
appropriately,  then: 


‘11. 


m-  + 


V  cos  a 


42. 


=  -  [cos  (qj,  n)(xi  mj)  +  sin  (qi,  „)(zi  m.)] 


(37a) 
(3  7b) 


(RHSF).  =  (x.  a 


1  P 


+  qi  «p  - 


cos  Cp  “  x^  sin  Gpinij 


2n 


[cos{q.^^)(zi  nij)  -  sm(qi^^)(xj  mj)](qi^j^  +  Gp)  -  Ki-  q- 


- 


4  «p 


V  cos  G, 


z  z  V 

O  O  00 

V~  ,,2 

00  V 


(x.  ip  -  q.)V^  (Xj  Gp  -  4i)(ip  sinGp) 

v2  V  2 

V  cos  V  cos  a 

00  p  00  p 

^^Gj  ^qj^^qi+n  “p^  RNi(Grp)  +  QDFj 


(37c) 
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M- 


V  cos  a. 


-  cos  (q;+n^(xj  m.)  sin  (q;+n)(z.  mj) 


(37d) 


"'22. 

1 


(37e) 


«P  +  «p 

2^ 


(RHSM)j  =  [(5^  _  xj  Op  -  q.  hp  cos  (qj^J  +  {x^  +  2q.  a 

-  xj  ap  sin  (qi+„)]{Xi  nij)  +  [(z^  ^  x.  -  qj  sin  (q^  +  n) 

■^)  cos  (qi_^,n)](Zi  mj)  -  " 

] 

(Xi  i  -  q.)V 


’'i  «p 


-  +  ^qjip  +  qi  «p  -  «p)  cos  (qj^^)](zj  mj)  -  .  qj 

-  AFj[x;(qj^J  -  (q^  -  q.)]  + 

(xi  ip  -  ij)(ip  sin  Hp) 


V 


V  cos  (z 


Y  cos 


z  z  V 

O  O  oo 

Y~  V  2 

OO  V 


+  (M-  +  ) 

^1  Hi 


X  (qi+n  +  ip)  +  PM;  (a.p)  +  QDMj 


(37f) 


Note  that  FN^  (o^p)  and  PMi(o'p')  are  the  "static"  aerodynamic  force  and 
moments,  respectively,  which  may  be  nonlinear  functions  of  and 
other  variables  as  well.  In  general,  curve  fitting  or  table  look-up 
interpolation  of  these  quantities  will  be  required. 

2.5.6  Numerical  Solution 


Equation  (35)  is  a  system  of  2n  second-order  differential  equations 
which  must  be  solved  simultaneously.  The  general  complexity  is  such 
that  a  numerical  solution  is  required.  To  accomplish  this,  a  computer 
program  (RKAM)  documented  in  Ref.  21  was  employed.  In  order  to 
apply  this  computer  program  written  for  simultaneous  first-order  dif¬ 
ferential  equations,  it  is  first  necessary  to  solve  Eq.  (35)  for  qj^  and 
qj^q_j^  explicitly.  This  yields 


"'ll.  "'22.  “  "'12.  "’21. 
11  11 


(38  a) 
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(38b) 


(39  a) 


(39b) 


Equation  (3  9)  can  be  reduced  to  a  system  of  4n  first-order  equations 
by  the  following  substitution; 

=  Sj^  k  =  1  —  2n  (40  a) 

Then 


S|^  =  k  =  1  —  2n 


(40b) 


Equation  (40)  is  a  set  of  4n  first-order  equations  to  be  solved  simul¬ 
taneously  subject  to  the  following  initial  conditions  from  Eq.  (37b); 


=  ®k  =  l---2n 

=  ^k  k  =  1  —  2n 


(40c) 


Equations  (40a),  (40b),  and  (40c)  are  in  the  form  required  for 
solution  by  the  RKAM  program  (Ref.  21).  The  resulting  computer  pro¬ 
gram  is  hereafter  referred  to  as  SMD2NDOF. 


The  sting  rigid  body  motion  variables  (up,  Up,  Up,  Xq,  Xq,  Zq,  and 
Zq)  must  be  known  as  a  function  of  time.  These  variables  are  required 
if  the  sting  system  moves  as  a  rigid  body  (for  example,  during^model 
to  tunnel  injection  or  sector  rotation).  and,  occasionally,  and 
q^o  must  be  known  as  a  function  of  time.  These  transient  variables  are 
currently  curve  fitted  by  a  Fourier  Series  or  a  least- squares  polynomial 
and  can  thus  be  computed  in  SMD2NDOF  when  needed  as  a  function  of 
time.  A  table  look-up  interpolation  representation  may  also  be  employed. 
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3.0  APPLICATIONS 


3.1  HIRT  STING  DESIGN  AND  TRANSIENT  RESPONSE 

References  22  and  23  provide  information  concerning  the  Advanced 
Transonic  Transport  or  the  Advanced  Technology  Transport  (ATT) 
model  which  is  a  representative  candidate  for  analysis.  Figure  8  illus¬ 
trates  the  predicted  history  of  the  dynamic  pressure  in  a  pilot,  high 
Reynolds  number,  Ludwieg-type  wind  tunnel.  This  information  is 
based  primarily  on  experience  with  such  a  VKF  pilot  facility  (Ref.  1). 
Maximum  steady-state  dynamic  pressures  (q^  )  were  selected  from 

conditions  reported  in  Ref.  22. 


Starting 

Time 


Figure  8.  Predicted  variation  of  dynamic  pressure  ratio  in  a  Pilot 
Ludweig-Tube  Wind  Tunnel. 


Both  the  steel  sting  and  the  C-S  sting  were  investigated  at  critical 
sting  design  load  conditions.  One  critical  point  was  test  condition  No.  7 
of  Ref,  22,  Table  2.  For  this  condition,  the  ATT  model  at  aq;  of  8.  0  deg 
generates  an  aerodynamic  normal  force  of  24,000  lb  (see  Fig.  3).  This 
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means  that  o-p  (the  initial  fixed  preset  angle  of  the  sting,  Fig.  6a)  plus 
(the  angular  deflection  at  the  model  location)  must  essentially  equal 
8  deg.  Thus,  in  order  to  initially  set  «p  before  the  run,  0i(t)  must  be 
computed. 

The  steel  sting  exhibited  highly  unsatisfactory  deflection  divergence 
characteristics  under  dynamic  load  conditions.  Just  a  small  initial 
preset  angle  of  attack  (up)  was  enough  to  cause  very  large  sting  deflec¬ 
tions  before  the  sting  strain  energy  restoring  force  became  sufficient 
to  prevent  further  deflection  (see  Fig.  9).  The  situation  is  analogous 
to  the  torsional  divergence  of  a  wing  (Ref.  11,  Chapter  3).  A  small 
increase  in  o-p  from  0.  60  to  about  0.  85  would  surely  result  in  catastro- 
pic  sting  failure.  The  situation  is  so  unstable  that  this  ATT  model 
could  not  be  tested  at  this  condition  with  a  steel  sting. 


Time,  t,  sec 


a.  ttp  =  0.60  deg 

Figure  9.  Motion  history  for  ATT  model  and  steel  sting,  condition 
No,  7,  basic  variation. 
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Time,  t,  sec 


b.  ap  =  0.25,  0.50,  and  0.75  deg 
Figure  9.  Concluded. 


Figure  10  illustrates  the  motion  history  {Bi  and  versus  time)  for 
the  ATT  model  and  the  composite  C-S  sting.  Here,  Up  is  5.0  deg,  and 
the  ATT  test  condition  No,  7  is  being  simulated.  Deflection  is  less  than 
half  of  that  observed  for  the  steel  sting  (Fig,  9a);  however,  this  is  a 
critical  condition  for  the  carbide  portion  of  the  sting  (x  <  48.  5  in. ) 
because  of  bending  stresses  at  these  loads.  Results  for  two  model 
weights  and  two  different  lumped-mass  representations  (n  =  1  and  n  =  15) 
are  shown  in  Fig,  10.  For  n  =  1,  note  that  essentially  the  same  results 
are  obtained  for  both  model  weights  (except  for  the  frequency).  The  fre¬ 
quencies  agree  very  well  with  predictions  for  the  C-S  sting  in  Fig.  4, 
particularly  the  MATVEC  results.  For  n  =  15,  the  sting  is  represented 
by  14  lumped  masses  so  that  sting  inertia  is  taken  into  account.  This 
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more  complex  simulation  of  the  sting  makes  about  15-percent  difference 
in  the  results  which  may  be  enough  to  warrant  inclusion  in  some  detailed 
analysis.  However,  the  n  =  15  results  shown  required  about  7  hr  of 
computation  time  on  the  IBM  370/165,  whereas  the  n  =  1  results  required 
only  a  few  minutes. 

Figure  11  illustrates  what  would  happen  if  the  dynamic  pressure 
were  a  step  input  for  ATT  test  condition  No.  7.  Results  for  two  model 
weights  are  shown  (1^.  was  held  constant).  Considerable  oscillation  is 
present  with  the  heavier  model  exhibiting  the  larger  magnitude.  It  is 
speculated  that  the  C-S  sting  would  fail  catastrophically  for  this  particu¬ 
lar  step  load  input  condition. 

Several  computer  runs  (not  illustrated  herein)  were  made  with  the 
basic  wind  tunnel  q^,  variation  (Fig.  8)  perturbed  so  as  to  yield  slightly 
faster  initial  rise  times.  The  resulting  sting  oscillation  amplitudes 
were  then  somewhat  more  than  for  the  basic  q^  variation  (Fig.  10)  but 
considerably  less  than  for  the  step  q^  input  (Fig.  11).  Both  the  heavy 
and  light  models  exhibited  practically  identical  oscillation  amplitudes 
(very  low)  for  the  basic  q^  variation  (Fig.  10,  n  =  1).  The  basic  q^ 

(Fig.  8)  aerodynamic  loading  is  then  essentially  like  a  "static"  loading 
since  inertial  loads  were  almost  nil.  For  the  step  input  (Fig.  11), 
the  loading  was  indeed  "dynamic"  and  the  heavy  and  light  models  did 
not  oscillate  at  the  same  magnitudes  because  of  the  strong  inertial 
effects.  It  is  worth  noting  that,  if  the  tunnel  dynamic  pressure  rise 
time  is  slower  than  the  basic  variation  of  Fig.  8,  then  the  steady- 
state  run  time  is  decreased  accordingly,  and  very  little  decrease  in 
sting  oscillation  amplitude  is  gained.  Consequently,  it  is  concluded 
that  the  example  wind  tunnel's  basic  q^  transient  variation  was 
essentially  optimum  from  the  viewpoint  of  both  sting  oscillation  and 
steady- state  run  time  for  model -sting  systems  of  this  size. 

In  lieu  of  actual  experimental  transient  data,  the  aerodynamic  loads 
for  the  ATT  were  generated  from  Eqs.  (22)  and  (26)  where  q^^,  =  qQo(t). 
Computer  runs  were  made  with  reasonable  estimates  of  C]^  ,  C]^ . ,  Cj^  , 
Cm  >  .  ,  and  Cj^  .  It  was  found,  for  this  particular  model  and  sting 

attachment  point,  th^  only  the  normal  force  (Cj^  )  influenced  the  sting 
deflection  to  any  practical  extent.  This  behavior  was  also  noted  in  a 
sting-model  motion  simulation  reported  in  Ref.  2.  These  stings  are  so 
stiff  that  aerodynamic  damping  exerts  practically  no  influence.  Grossly 
exaggerated  static  moment  inputs  (both  stable  and  unstable  via  Eq.  (28)) 
did  produce  noticeable  changes  in  deflection  (approximately  T  15  percent, 
respectively). 


38 


AEDC-TR-76-41 


Figure  10.  Motion  history  for  ATT  model  and  carbide-steel  sting, 
condition  No.  7,  basic  variation,  Up  =  5.0  deg. 

Model 

Sym  Weight,  lb  u ,  cps  Remarks  n  DOF 

-  350  10.0  CNo  Only  1  2 

-  1,270  6.0  Cnq  Only  1  2 


Figure  11.  Motion  history  for  ATT  model  and  carbide-steel  sting, 
condition  No.  7,  ap  =  5.0  deg. 
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3.2  VKF  TUNNEL  F  STING  DYNAMICS 


Recent  results  obtained  in  VKF  Tunnel  F  have  provided  an  opportunity 
for  comparison  with  experimental  motion  data.  During  this  test  program, 
the  sting  angular  deflection 'was  found  to  be  greater  than  could  be  account¬ 
ed  for  by  model  aerodynamic  loads  alone.  High-speed  motion  pictures 
(framing  rate  »4000  frames/sec)  were  taken  during  several  shots  to 
provide  a  record  of  the  sting-model  deflection.  It  was  speculated  that 
there  were  significant  aerodynamic  and  inertial  loads  on  the  sting  itself 
which  contributed  to  the  observed  deflection.  Consequently,  the  SMD2NDOF 
computer  program  was  applied  to  simulate  the  sting  motion,  both  with 
and  without  sting  aerodynamic  loads.  A  general  description  of  the  VKF 
Hypervelocity  Tunnel  F  is  contained  in  Ref.  24. 

The  Tunnel  F  sting-balance  was  modeled  structurally  via  the 
STING -1  program.  Some  experimental  calibration  rotational  deflection 
data  for  the  balance  alone  and  the  balance  plus  sting  were  required  to 
obtain  a  reasonable  stiffness  representation  of  this  model  support 
system.  The  basic  sting  utilized  for  this  particular  test  is  not  as  stiff 
as  its  exterior  dimensions  woiild  suggest  because  of  numerous  joints 
and  cross-sectional  area  discontinuities.  Figure  12a  illustrates  the 
actual  sting  used  during  the  tests,  including  the  force  balance.  The 
general  complexity  of  this  actual  sting  force  balance  system  is  such  that 
a  number  of  dimensional  iterations  were  necessary  via  the  STING -1  pro¬ 
gram  to  obtain  the  structural  stiffness  synthesis  illustrated  in  Fig.  12b. 

A  somewhat  more  complete  synthesis  of  the  sting-force  balance  system 
is  feasible  and  should  be  incorporated  in  future  analyses. 

Transient  model  loads  (corrected  for  inertia  effects)  were  measured 
during  the  test.  These  were  curve  fitted  with  a  Fourier  series  and 
input  as  forcing  functions  (N'pj(t)  and  Mq72(t))  to  the  SMD2NDOF  program 
(Fig.  13).  Strictly  speaking,  the  forward  and  aft  normal  loads  measured 
by  the  force  transducers  in  the  balance  were  converted  into  a  normal 
force  and  pitching  moment  which  were  assumed  to  act  at  a  single  model- 
balance  attachment  point  (the  physical  eg  location  of  the  model.  Fig.  12b). 
This  pitching  moment  makes  a  sizable  contribution  to  the  sting  deflection 
and  cannot  be  neglected. 

Figure  14  illustrates  the  Tunnel  F  transient  flow  velocity  history 
for  this  particular  test. 
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Figure  13.  Transient  aerodynamic  loads  on  the  Tunnel  F  model. 


Figure  14.  Transient  flow  velocity  history  for  VKF  Tunnel  F. 
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Since  airloads  on  the  sting  were  to  be  included,  the  sting  was  sub¬ 
divided  into  14  sections  (concentrated  masses,  n  =  15,  including  the 
model).  No  experimental  data  were  available  for  the  aerodynamic  loads 
acting  directly  on  the  sting  itself.  Consequently,  these  airloads  were 
estimated  in  the  following  manner.  It  was  assumed  that  the  airloads 
(NTj^)  on  the  individual  sting  masses  (mj^)  had  the  same  timewise  varia¬ 
tion  as  the  airload  (Nrp^)  on  the  model.  However,  the  magnitude  of  the 
sting  airloads  is  not  equal  to  model  airloads  because  of  different  plan- 
form  areas  and  different  shape  effects.  The  sting  is  composed  of  a 
series  of  circular  cylinders  and  cone -frustums,  while  the  model  is 
essentially  conical.  The  airload  normal  force  on  a  cylindrical  concen¬ 
trated  mass  is  then 


{41a) 


where  (model)  and  C]\|-^  (sting  m^)  are  based  on  planform  area. 
for  the  sting  was  estimated  from  the  Newtonian  normal  force  coefficient 
for  circular  cylinders  (Ref.  25,  Chapter  17).  The  ratio  (Cjg./C];^^) 
remains  practically  constant  from  15  to  17  deg  which  essentially  brackets 
the  observed  angles  of  attack.  A 
the  model  and  m^,  respectively. 

the  actual  sting  (Fig.  12a)  were  employed.  For  cone-frustum  sections 
of  the  sting,  shape  effects  were  assumed  to  be  similar  to  those  of  the 
model  so  that 


and  Ap.  are  the  planform  areas  of 
Kanform  areas  (Ap J  corresponding  to 


'  \  Pi/ 


(41b) 


Equations  (41a)  and  (41b)  have  the  general  form: 


Njit) 


C,  Nt  (t) 


(41c) 


The  loads  on  the  sting  section  were  assumed  to  act  at  the  lumped  mass 
eg  and  produce  no  pitching  moment  (Mq;'^(t)). 

Only  the  last  12  sting  masses  were  considered  to  be  aerodynami- 
cally  loaded  since  the  two  forward  portions  of  the  sting  are  inside  the 
model.  The  summation  of  C^  (i  =  4  to  15)  yields  the  ratio  of  the  sting 
total  normal  airload  to  the  model  normal  airload.  In  this  particular 
case, 

2  C.  =  3.31  (42) 

i=4  ' 


\ 
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This  indicates  that  the  sting  airloads  are  about  three  times  the  model 
airloads.  Fortunately,  the  sting  loads  are  distributed  primarily  along 
the  stiffen  aft  sections  which  alleviates  their  deflection  contribution 
considerably.  The  theoretical  magnitude  of  the  sting  airload  contri¬ 
bution  to  the  static  deflection  can  be  computed  via  the  deflection  influ¬ 
ence  coefficient  approach  (Section  2.3).  Theoretical  results  from 
SMD2NDOF  indicate  that 

01  (with  sting  loads) 

o  -7-^-^ - - -s  «  f-12 

0j  (without  sting  loads) 

Figures  15  and  16  illustrate  both  experimental  results  and 
SMD2NDOF  predictions  for  this  Tunnel  F  test.  Motion  predictions 
were  made  with  and  without  sting  aerodynamic  loads.  Each  prediction 
requires  about  45  min  actual  computation  time  on  the  IBM  3  70/165. 

The  experimental  data  (from  high-speed  motion  pictures)  are  9-^  (Fig. 

15)  and  z^ose  ^6)  versus  time.  Zj^Qgg  was  computed  by  SMD2NDOF 

as  illustrated  in  Fig.  17. 

The  predicted  dynamic  oscillation  amplitude  (oscillation  about  the 
mean  or  static  level.  Fig.  15a)  is  somewhat  greater  than  the  experi¬ 
mental  magnitude,  particularly  near  the  end  of  the  run.  It  is  evident 
from  the  experimental  data  (Figs.  15a  and  16a)  that  considerable  damp¬ 
ing  of  the  oscillation  is  occurring.  From  previous  results  (Section  3.  1), 
it  is  known  that  the  damping  is  structural  rather  than  aerodynamic. 

Sting  structural  damping  was  included  in  the  SMD2NDOF  predictions 
illustrated  in  Figs.  15b,  d  and  16b,  d.  The  numerous  joints  of  the  actual 
sting  (Fig.  12a)  make  an  important  contribution  to  structural  damping 
(Section  2.  5.  4.  2).  This  is  fortunate  from  a  "static"  data  accuracy 
viewpoint  and  more  than  compensates  for  the  apparent  decrease  in 
stiffness  introduced  by  the  joints. 

The  average  magnitude  (mean  value  of  the  amplitude)  of  B-y  versus 
time  is  shown  in  Fig.  15c.  Examination  of  these  mean  amplitudes 
reveal  a  maximum  difference  of  about  0.  08  deg  between  experiment 
and  prediction  with  sting  loads,  and  about  0.  20  deg  difference  between 
experiment  and  prediction  without  sting  loads.  The  predicted  results 
with  sting  loads  included  are  in  remarkably  good  agreement  with  ob¬ 
served  magnitudes  near  the  beginning  of  the  run.  These  results  show 
that  it  is  feasible  to  account  for  sting  aerodynamic  loads  on  0^  (and 
hence 

There  is  a  slight  discrepancy  between  the  predicted  (48  cps)  and 
experimental  (55  cps)  oscillation  frequencies.  When  modeling  the 
sting  balance  system,  the  structural  stiffness  synthesis  was  of  primary 
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concern,  and  when  those  criteria  were  satisfied,  coi  was  48  cps  for 
sting-balance  plus  model  (STING- 1  program).  This  was  deemed  satis¬ 
factory  since  the  system  frequency  was  known  to  be  on  the  order  of  50 
cps.  A  large  amount  of  structural  damping  distorts  the  sinusoidal 
amplitude  variation  considerably  (Figs.  15b  and  16b)  and  some  discre¬ 
tion  is  required  when  computing  the  frequency  from  a  limited  number 
of  cycles  of  experimental  data. 

Sym 


a.  Transient  amplitude  without  structural  damping 


b.  Transient  amplitude  with  structural  damping 
Figure  15.  Theoretical  and  experimental  angular  deflections  of  the  Tunnel  F  model. 
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Sym 

o  Experimental  High-Speed 


d.  Oscillation  andplitude 
Fi^re  15.  Concluded. 


46 


AEDC-TR-76-41 


Arc  Fired 


-  S(V\D2ND0F  (n  =  15)  with  Sting  Airloads, 

No  Structural  Damping 

- SMD2ND0F  (n  =  15)  without  Sting 

Airloads,  No  Structural  Damping 

o  Experimental,  High-Speed  Motion-Picture 
Data,  Q  ~  55  cps 


Time,  sec 

a.  Transient  amplitude  without  structural  damping 

Arc  Fired 

A  ^  ^ 

A  A 

SMD2ND0F  with  Sting  Airloads 
and  Structural  Damping, 

FCOFj  =  FCOMj  =  0. 050 

Experimental  High-Speed 
Motion-Picture  Data 

Time,  sec 


Figure  16.  Theoretical  and  experimental  translational  model  nose  tip  deflections. 
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c.  Average  amplitude 


o  Experimental  High-Speed 


d.  Oscillation  amplitude 
Figure  16.  Concluded. 
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A  realistic  assesment  of  the  magnitude  in  Tunnel  F  is  of 

more  than  academic  interest  since  it  is  difficult  to  measure  experimen¬ 
tally  and  it  does  contribute  to  the  total  transient  angle  of  attack  of  the 
model  (urp^).  From  Eq.  (29), 


=  «p  + 


+  a 


plunge 


(43) 


where  in  the  present  case, 

®plunge^  =  ~  Y  cos  a  (44) 

oo  p  '  ' 

Theoretical  results  (including  sting  aerodynamic  loads  and  struc¬ 
tural  damping)  for  both  Q'piunge  illustrated  in  Fig.  18.  The 

maximum  value  of  the  predicted  u'piunge  approximately  one  tenth  of 
a  degree.  Because  of  the  large  amount  of  structural  damping  present, 
it  is  concluded  that  the  maximum  ^pi-Q^ge  will  normally  be  less  than 
0.  025  deg  for  force  tests  in  Tunnel  F  when  steady-state  data  are  taken 
after  the  initial  starting  transients. 


49 


Sym 


Time,  sec 

Figure  18.  Theoretical  results  for  Opiunge  and  model  eg  translation  {zi ). 

4.0  SUMMARY 

This  investigation  of  wind  tunnel  sting  structural  dynamics  has 
essentially  followed  two  paths.  These  are: 


4.1  STING  STRUCTURAL  DESIGN 

A  brief  analysis  (Section  2.  1)  has  shown  that,  for  a  given  wind 
tunnel  facility  (with  fixed  dimensions),  any  real  gain  in  sting  stiffness 
must  ultimately  come  from  inherent  stiffness  properties  (E)  of  the  sting 
material  itself.  This  means  that  materials  stiff er  than  steel  must  be 
considered.  Certain  carbide  alloys  (sintered  tungsten -titanium 
carbide)  are  about  three  times  stiffer  than  steel  and  warrant  serious 
attention  as  a  sting  construction  material.  However,  any  sting  material 
must  be  able  to  withstand  bending  stresses  of  about  120,  000  Ib/in. 

It  is  not  clear  that  carbide  alloys  in  sting  size  beams  will  consis¬ 
tently  handle  bending  stresses  of  this  magnitude.  This  should  be  checked 
by  standard  experiments. 


AEDC-TR-76-41 


To  assist  in  structurally  analyzing  realistic  sting  configurations, 
a  computer  program  (STING- 1)  was  written  (Section  2.2).  With  geome¬ 
trical  and  material  properties  as  input  data,  STING- 1  computes  deflec¬ 
tions,  first  natural  frequency,  and  influence  coefficients  (spring  con¬ 
stants)  for  cantilever  beams  with  mixed  tapered  and  untapered  sections. 
The  material  of  each  section  is  arbitrary  so  that  a  certain  type  of  com¬ 
posite  sting  may  be  analyzed.  Sting  design  applications  are  obvious; 
however,  perhaps  the  most  important  application  so  far  has  been  the 
computation  of  structural  influence  coefficients  (Section  2.3)  required 
for  certain  dynamic  analysis  procedures  (Section  2.4  and  2.  5.3).  The 
STING- 1  program  is  based  primarily  on  elementary  strength  of  material 
concepts  but  is  general  and  flexible  enough  to  model  most  of  the  more 
common  sting-balance  configurations. 

The  MATVEC  computer  program  (Ref.  13)  which  solves  the  standard 
matrix  eigenvalue  problem  is  recommended  for  the  computation  of  sting 
steady-state  vibration  natural  frequencies  and  mode  shapes.  The  input 
required  is  the  [CTF^j  m^]  matrix  (Section  2.4),  where  both  CTF^j 
and  m.  may  be  computed  by  the  STING-1  program.  Prior  experience 
in  a  structural  laboratory  with  cantilever  beams  has  verified  MATVEC 
predictions. 


4.2  STING  STRUCTURAL  DYNAMICS 


Lagrangian  mechanics  has  been  employed  to  derive  the  coupled  dif¬ 
ferential  equations  of  planar  motion  for  a  rather  general  sting-balance 
system  (Section  2.  5).  This  system  is  represented  by  concentrated 
masses  which  are  tied  together  by  spring  constants  (influence  coefficients) 
and  may  be  subjected  to  arbitrary  known  transient  forces  and  moments. 
The  formulation  allows  sting  loads  as  well  as  model  loads  to  be  included. 
Also  the  system  is  allowed  to  move  (rotate  and  translate)  as  a  rigid  body. 
Thus  with  the  proper  input,  transient  conditions  occurring  during  model 
to  tunnel  injection  and/or  sector  rotation  may  be  investigated. 

The  general  formiilation  is  such  that  a  numerical  solution  is  required 
for  all  but  the  simplest,  almost  trivial,  cases.  Consequently,  a  stand¬ 
ard  program  (Ref.  21)  for  solving  simultaneous  differential  equations 
was  employed  and  modified  as  required.  The  resulting  program 
(SMD2NDOF)  is  versatile  enough  to  simulate  a  wide  variety  of  sting 
structural  dynamic  situations.  The  sting-model  synthesis  does  not  have 
to  be  an  elaborate  multi-mass  system  to  provide  useful  information.  For 
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example,  see  Section  3.  1  where  a  single  mass  representation  was  em¬ 
ployed  for  the  analyzed  sting-model.  The  simulation  of  the  Tunnel  F 
system  (Section  3.  2)  required  a  more  elaborate  multi-mass  representa¬ 
tion  since  sting  aerodynamic  loads  were  to  be  included. 
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NOMENCLATURE 


Three  basic  axis  systems  are  used  in  the  sting  planar  motion  anal¬ 
ysis  of  Section  2.  5.  Figure  6  illustrates  these  axis  systems,  and 
they  are  further  explained  as  follows: 


Xjj,  Zjj 


This  is  a  fixed  inertial  axis  with  Xj  aligned 
parallel  to  The  location  of  the  origin  is 

arbitrary. 

This  axis  is  fixed  to  the  undeflected  sting  and 
moves  (rotates  and  translates)  with  the  unde¬ 
flected  sting.  The  origin  is  at  the  center  of 
the  sector  rotation  and  Xjj  is  directed  along 
the  undeflected  sting  centerline.  This  is  the 
basic  axis  system  from  which  sting  deflections 
are  measured. 

This  axis  system  is  fixed  to  the  i^^  lumped  mass 
of  the  system  and  rotates  and  translates  with 
this  mass,  tangent  to  the  deflected 

sting  centerline.  The  origin  is  arbitrary  but 
if  chosen  at  the  mass  eg  will  eliminate  terms 
involving  %m^  and  in  the  motion  equations. 


Additional  symbols  and  nomenclature  and  abbreviations  used  herein 
are  defined  as  follows: 


AF^ 

Api 

^Pi 

ATT 


ak 

t>k 

CAi 


Instantaneous  axial  force  on  the  model,  lb, 

(Eq.  (32)  and  Fig.  8) 

Planform  area  of  model,  ft^ 

Planform  area  of  the  i^^  mass  of  the  sting,  ft^ 

Advanced  transonic  transport  or  advanced  technol¬ 
ogy  transport. 

Initial  value  of  qk  at  t  =  tQ  (Eq.  (39b)) 

Initial  value  of  qk  at  t  =  tQ  (Eq.  (3  9b)) 

Instantaneous  axial-force  coefficient  (Eq.  (32)) 

Angular  deflection  at  i  caused  by  unit  force  at 
j  (Eq.  (7)) 
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CAM^j 

Ci 

Cm 

Cm 

mq 

Cm'p 

Cm 

mQ, 

C 

Cn 

CNq 

Cnt 

Cn^ 

Cn. 

a 

C-S 

CTF.. 

13 

CTM.. 

13 


eg 


Angular  deflection  at  i  caused  by  unit  moment  at 
3  (Eq.  (7)) 

Constant  defined  by  Eq.  (41) 

Flexibility  influence  coefficients  (Eq.  (6)) 

Static  aerodynamic  pitching -moment  coefficient 
(Eq.  (28)) 

Aerodynamic  pitch  damping  moment  coefficient 
caused  by  q^^  (Eq.  (2  7)) 

Total  aerodynamic  pitching -moment  coefficient 
(Eq.  (27)) 

Slope  of  static  pitching-moment  coefficient, 
dC^/dQ-  (Eq.  (28)) 

Aerodynamic  pitch  damping  moment  coefficient 
caused  by  a  (Eq.  (27)) 

Static  aerodynamic  normal-force  coefficient 
(Eq.  (24)) 

Aerodynamic  normal-force  coefficient  caused  by 
qj^  (Eq.  (23)) 

Total  aerodynamic  normal-force  coefficient 
(Eq.  (23)) 

Slope  of  static  normal  force  coefficient, 
dC^lda  (Eq.  (24)) 

Aerodynamic  normal  force  coefficient  caused 
by  a  (Eq.  (23)) 

Carbide -Steel  composite  sting 

Translational  deflection  at  i  caused  by  unit 
force  at  j  (Eq.  (7)) 

Translational  deflection  at  i  caused  by  unit 
moment  at  j  (Eq.  (7)) 

Distance  from  beam  neutral  axis  to  surface, 

(Eq.  (2)),  c  =  d/2  for  circular  cross-section 
stings 

Center  of  gravity  or,  more  correctly,  center  of 
mass  location 


DOF 


Degrees  of  freedom  (number  of  independent 
coordinates  necessary  to  describe  motion) 
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d 

E 

El  to  E5 

F 

FN 

FCOF^, 

FCOMj^ 

fk 

g 

I 

ID 

V 


i, 

K 

K.. 

ij 

L 

Ll  to  L5 
Nog 
M 


M 


T 


M 


M- 

a 


Aerodynamic  reference  length,  ft  or  sting 
diameter,  in. 

Modulus  of  elasticity  or  Young's  Modulus, 
lb/ in.  ^ 

E  for  beam  sections  in  STING- 1  computer 
program  (Fig.  2) 

Friction  force  (Eq.  (33)) 

Static  aerodynamic  normal  force,  lb 
(FN  =  S  C^) 

Structural  damping  constants  (Eq,  (34)) 

General  functional  form  of  (Eq.  (39a)) 

Gravitational  constant,  Ibm/lbf /ft/sec^ 

Area  moment  of  inertia  of  beam  cross  section, 
in.  4  or  ft‘^ 

Internal  diameter  of  a  hole  through  a  sting 
section,  in.  (Fig.  12) 

Mass  moment  of  inertia  of  the  model  or  a 
discrete  lumped  mass  of  the  system  represen¬ 
tation,  ft-lb-sec2  (Eq.  (17d)) 

General  indices 

Linear  spring  constant,  force /deflection 

Stiffness  influence  coefficients  (Eq.  (8a)) 

Length  of  sting,  in.  or  ft 

Length  of  individual  sting  sections  (Fig.  2) 

Length  from  model  nose  tip  to  eg  position,  in. 
(Figs.  12  and  17) 

Total  bending  moment  acting  on  sting  at 
station  x,  in.  -lb,  or  ft -lb 

Total  aerodynamic  pitching  moment,  ft-lb 
(Eq.  (26)) 

(q„S  d2/2  VJ 

q 

(q„S  dh2  VJ  C 

III  • 

a 
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m 

Mass,  lb-sec2/ft 

rrii 

Mass  of  system  component,  lb-sec^ /ft 

"^lli 

Quantity  defined  by  Eqs,  (36a)  or  (37a) 

nil2i 

Quantity  defined  by  Eqs.  (36b)  or  (37b) 

^21i 

Quantity  defined  by  Eqs.  (36d)  or  (37d) 

m22i 

Quantity  defined  by  Eqs.  (36e)  or  (37e) 

N 

Number  of  DOF,  N  =  2n  in  the  analysis  of 

Section  2.  5 

Total  aerodynamic  normal  force,  lb  (Eq.  (22)) 

N1  to  NEND 

Station  numbers  in  STING- 1  program  (Fig.  2) 

Nq 

(q.,  S  d/2  V^) 

tn 

{q^  S  d/2  V„) 

n 

Number  of  lumped  or  concentrated  masses 
employed  in  system  representation 

OD 

Outer  diameter  of  sting,  in. 

PM 

Static  aerodynamic  pitching  moment,  in.  -lb  or 
ft-lb  (PM  =  q^  Sd  Cm) 

Qk 

Non- conservative  generalized  force  (or  moment) 
(Section  2.5.4) 

QDP. 

General  structural  damping  force,  lb  (Eq.  (34a)) 

QDMi 

General  structural  damping  moment,  ft-lb 
(Eq.  (34b)) 

Generalized  independent  deflection  coordinate, 
used  also  with  various  indices  such  as  i,  j,  or 
i  +  n 

% 

dqj^/dt 

d^q^/dt^ 

Qr 

Up  +  9,  rotational  portion  of  ur^,  radian/ sec 

Qoo 

Free -stream  dynamic  pressure  which  may  be  a 
transient  quantity,  Ib/ft^ 

Maximum  value  of  q^^,  Ib/ft^ 

RHl  to  RH5 

Hole  radius  in  sting  sections,  STING- 1  Program 
in.  (Fig.  2) 
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(RHSF). 

Right-hand  side  of  i^^  force  Eq.  (35a) 

(Eqs.  (36c)  or  (37c)) 

(RHSM)^ 

Right-hand  side  of  i^^  moment  Eq.  (35b) 

(Eqs.  (36f)  or  (37f)) 

R(N1)-R{NEND) 

Radius  at  sting  section  discontinuities  in 
STING- 1  Program,  in.  (Fig.  2) 

S 

Reference  area,  ft^ 

Sk 

qj^,  (Eq.  (40a)) 

Sk 

%,  (Eq.  (40b)) 

T 

Total  kinetic  energy  of  sting -model  system, 
ft -lb  (Section  2,5.2) 

T. 

1 

th 

Kinetic  energy  of  the  i  mass  of  the  system 
representation,  ft -lb  (Section  2.  5.  2) 

t 

Time,  sec 

Initial  time,  usually  taken  as  zero  before 
transient  conditions  begin,  sec 

U 

Total  potential  energy  of  the  sting-balance 
system,  ft -lb  (Section  2.  5.  3) 

^ABS 

Absolute  velocity  of  a  particle  of  mass  (dm), 
ft/ sec  (Section  2.5.2) 

V., 

Free-stream  tunnel  veolcity  which  is  allowed 
to  be  a  transient  quantity,  ft/sec.  must 

be  known  to  compute  Q'piunge- 

Wi 

Weight  of  a  concentrated  mass,  lb 

Wm 

Weight  of  model,  lb 

Wl  to  W5 

Specific  weights  of  sting  section  material  in 
STING- 1  program,  Ib/in.  ^  (Fig.  2) 

w 

Specific  weight  of  sting  material,  lb/ in.  ^ 

X 

General  coordinate  along  undeflected  sting 
centerline  measured  from  sector  rotation 
point,  in. ;  also  general  distance  coordinate 

^0 

Horizontal  location  of  Zjj  axis  origin  with 

respect  to  the  Xj,  Zj  axis  origin,  ft 
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X 


O 


X,- 


X 


Hi 


xill- 


z 


z 


o 


o 


z 


i 


Horizontal  velocity  component  of  Xjj,  Zjj 
axis  (sting-model  system  moving  as  a  rigid 
body),  ft/ sec 

Horizontal  acceleration  component  of  Xjj,  Zjj 
axis,  ft/sec^ 

xjj.,  location  of  i'th  mass  measured  along  unde¬ 
flected  sting  centerline  from  model  attachment 
point,  in.  (Eq.  (16c),  Fig.  6) 

Location  of  eg  of  mass  (mj^)  with  respect  to 
Xm.  ,  Ziii.  origin,  ft  (Eq.  (17b),  Fig.  6c) 

Horizontal  location  of  particle  of  mass  (dm)  with 
respect  to  Xj,  Zj  axis  origin,  ft  (Eq.  (14a), 

Fig,  6c) 

Xj^,  location  of  mass  measured  along  unde¬ 
flected  sting  centerline  along  Xu  from  Xn, 

Zjj  axis  origin,  ft  (Eq.  (16c),  Fig.  6) 

Location  of  particle  of  mass  (dm)  with  respect 
to  Xjjj. ,  origin,  ft,  (Fig.  6b) 

General  deflection  coordinate  of  sting  center- 
line,  measured  perpendicular  to  undeflected 
sting  centerline,  in. 

Vertical  location  of  Xjj,  Zjj  axis  origin  with 
respect  to  the  Xj ,  Zj  axis  origin,  ft 

Vertical  velocity  Xjj ,  Zjj  axis  (sting-model 
system  moving  as  a  rigid  body),  ft/sec 

Vertical  acceleration  component  of  Xjj,  Zjj 
axis,  ft/sec^ 

tb 

Zji  or  qj^,  deflection  of  i  mass  relative  to 
undeflected  sting  centerline,  in.  (Eq.  (16a), 

Fig.  6)  also  normalized  deflections  (eigenvectors) 
in  Eq.  (11) 

dzj_/dt  or  ,  ft/sec 

d^Zj^/dt^  or  qj^,  ft/sec 

Location  of  eg  of  mass  (mj^)  with  respect  to 
XjTT^ ,  Zjjj^  axis  origin,  ft  (Eq.  (17c),  Fig.  6c) 


60 


AEDC-TR-76-41 


2^nose 


^amplitude 


^plunge 


Vertical  location  of  particle  of  mass  (dm)  with 
respect  to  Xj,  Zj  axis  origin,  ft  (Eq,  (14b), 

Fig.  6c) 

Zj^  or  ,  deflection  of  i^h  mass  with  respect 
to  Xjj  axis,  in,  (Eq.  (16a),  Fig.  6) 

Location  of  particle  mass  (dm)  with  respect 
to  ,  Zjjj^  axis  origin,  ft  (Fig.  6b) 

Average  amplitude  of  z,  ft  or  in. 

Maximum  amplitude  of  Zj_ ,  ft  or  in. 

Amplitude  of  model  nose  tip  translation,  in. 

(Fig.  17) 

Oscillation  amplitude  of  z,  ft  or  in.  , 

^^amplitude  ~  ^max  ~  ^av^ 

Pitch  angle  of  sting-model  acting  as  a  rigid 
system  (rotation  of  X^j ,  Zjj  axis),  deg  or  rad 

Angular  rate  of  Xjj,  Zjj  axis  rotation,  rad/sec 

Angular  acceleration  of  Xjj,  Zjj  axis,  rad/sec^ 

Angle  of  attack  caused  by  transverse  plunging 
velocity  components,  deg  or  rad  (Eqs.  (29),  (43), 
and  (44)) 

Total  angle  of  attack,  deg  or  rad  (Eq.  (29)) 

dttrp/dt,  rad/sec  (Eq.  (30)) 

Moment  arm  of  axial  force,  ft  ( ~ 
zi  +  z^)  (Fig,  7) 

Maximum  bending  strain  at  station  x,  in/ in.  , 

(Eq.  (3)) 

9i+n »  rotational  deflection  of  sting  at  station  x^^, 
deg  or  rad  (Eq.  (16b)) 

^i+n’  ^ad/sec 

*^i+n'  ^ad/sec^ 

Average  amplitude  of  0^^,  deg 

Maximum  amplitude  of  ,  deg 
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9 


^amplitude 


A 


Oscillation  amplitude  of  9^,  deg 


imax 


-  01  ) 
■'•av 


1/w^,  (CTFjj  m^)  matrix  eigenvalue,  (Eq,  (11)) 

o 

Maximum  bending  stress  at  station  x,  Ib/in. 
(Eq.  (2)) 


Material  yield  stress,  Ib/in.  2  (Table  1) 

Frequency  of  vibration  of  the  i^^  mode,  cps 
(i  =  1,2,3,  etc. ) 
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